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I.  INTRODUCTION 


A.  Background 

This  study  is  conducted  within  the  framework  of  the  detection 
problem  encountered  in  an  active  sonar  employing  multiple  receiving  ele¬ 
ments.  Active  sonar  is  a  device  for  discovering  distant  objects  such  as 
submarines,  ships,  fish,  etc.  It  utilizes  a  transducer  which  generates 
acoustic  energy  of  wavelengths  on  the  order  of  1  to  20  cm  in  the  form  of 
pulses  of  large  amplitude  and  brief  duration.  When  this  energy  strikes 
a  reflecting  object  it  is  scattered  in  all  directions.  Some  of  this 
energy  returns  to  the  transducer  as  signals  which  are  received  by  the 
sonar  receiver.  This  study  considers  the  manner  in  which  these  signals 
are  processed  by  the  sonar  receiver. 

At  first  glance  one  might  think  that  targets  could  be  detected  at 
any  range  merely  by  increasing  the  amplification  of  the  received  signals. 
In  actual  sonar  conditions,  however,  there  are  several  types  of  inter¬ 
ference  which  are  detrimental  to  the  sonar  receiver's  performance.  These 
sources  of  interference  can  be  classified  as  system  noise,  ambient  noise, 
and  reverberation  noise. 

System  noise  includes  both  thermal  noise  produced  in  the  receiver 
itself  and  noise  produced  by  the  motion  of  the  array  through  the  water. 

In  this  study  flow  noise  is  not  significant  since  a  stationary  array  is 
used.  Also,  the  other  noise  sources  are  assumed  to  be  the  predominant 
effects  which  limit  the  receiver's  performance. 

Ambient  noise  sources  are  often  quite  powerful  and  contain  a  wide 
range  of  frequencies.  Some  examples  of  this  type  of  noise  source  are 
wave  action  at  the  surface,  meteorological  effects,  and  noise  generated 
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by  distant  ships.  Ambient  noise  sources  are  of  particular  interest  to 
the  designer  of  passive  sonars.  In  this  study  the  ambient  noise  field  is 
assumed  to  be  negligible  as  compared  to  reverberation  noise. 

Reverberation  is  the  noise  source  which  poses  the  most  difficult 
problem  for  an  active  sonar  receiver.  Since  this  noise  is  caused  by  the 
transmitted  signal,  it  has  many  characteristics  similar  to  the  echo  from 
an  underwater  object.  It  does  no  good  to  increase  the  gain  of  the 
receiver's  amplifier  because  the  power  of  the  reverberation  is  increased 
proportionately.  There  are  three  types  of  reverberation:  volume,  bottom, 
and  surface  reverberation,  and  these  are  caused  by  reflection  of  the 
transmitted  signal  from  irregularities  in  the  body  of  water  and  from  the 
lower  and  upper  surfaces  of  the  water  respectively.'*'  The  type  of  noise 
considered  in  this  study  is  surface  reverberation. 

In  active  sonar  signal  detection  the  target  echo  is  often  received 
at  a  low  level  compared  to  the  energy  in  the  surface  reverberation.  In 
an  attempt  to  enhance  the  target  signal,  some  current  Navy  sonars  employ 
matched  filtering  to  process  the  received  sonar  signal.  Since  this  type 
processing  consists  of  merely  correlating  the  incoming  signal  with  a 
replica  of  the  transmitted  waveform,  it  is  optimum  (in  the  Nevman-Pearson 
sense)  when  the  noise  is  a  Gaussian  process,  both  independent  and  sta¬ 
tionary.  However,  surface  reverberation  is  in  general  a  nonstationary 
process  with  a  finite  correlation  time;  therefore  a  processor  that  must 
cope  with  surface  reverberation  must  be  more  complicated.  There  is  an 
extensive  body  of  literature  on  tneoretica>  signal  processing  techniques 
designed  to  determine  whether  or  not  a  target  echo  is  present  in  a  given 
sonar  return.  We  have  some  knowledge  of  more  optimum  ways  to  process 
these  signals,  but  we  have  no  fee.  for  tne  efficiency  of  those  improved 
processing  techniques.  Both  a  brief  theoretical  review  and  some  experi¬ 
mental  effort  seem  justified.  Specifically,  it  would  be  valuable  to 
know  how  a  given  signal  processing  technique  performs  in  the  presence  of 
surface  reverberation. 


B. 


Previous  Studies 


As  stated  previously,  much  work  has  already  been  done  to  develop 

better  techniques  for  processing  sonar  signals  from  multiple  receiving 
2 

elements.  Mermoz  has  developed  the  optimum  receiver  structure  to  maxi¬ 
mize  the  output  signal-to-noise  ratio.  In  that  work,  the  noise  process 
is  assumed  stationary  with  a  known  correlation  matrix,  but  otherwise 
arbitrarily  distributed.  The  result  is  a  set  of  generalized  matched 
filters.  Each  filter  response  depends  on  the  signal  and  the  matrix  of 
spectral  densities  of  the  noise  process.  An  excellent  summary  of  Mermoz 's 
work  is  given  by  Horton.'*' 

3 

Middleton  and  Groginsky  specified  the  receiver  structure  which  is 

optimum  in  the  Bayes  sense,  i.e.,  minimum  average  cost  to  the  observer. 

That  paper  dealt  with  Gaussian  signals  in  the  presence  of  Gaussian  noise, 

but  stationarity  was  not  an  assumption.  One  of  the  more  important  results 

presented  was  that  factorization  of  the  spatial-temporal  operations  is 

not,  in  general,  possible  in  optimum  systems.  Factorization  here  is  the 

"separation  of  the  required  filter  into  two  successive  operations — the 

first  depending  only  on  the  geometry  of  the  array,  the  second  depending 

3 

only  on  the  statistics  of  the  noise  process."  Current  Navy  sonars  uti¬ 
lize  factored  processing  by  first  beamforming  and  then  performing  some 
temporal  operation  such  as  matched  filtering.  The  optimum  detector 
developed  in  that  paper  utilizes  a  set  of  time  varying  filters  which  are 
determined  by  the  covariance  functions  of  the  noise  and  signal  processes. 

In  both  of  these  works  the  covariance  functions  of  the  noise  process 
are  assumed  known.  For  surface  reverberation,  the  process  that  is  a 
major  limitation  of  active  sonar  receivers,  this  assumption  is  hardly  a 
minor  one.  In  recent  years,  however,  several  experimental  studies  have 

been  conducted  which  considered  the  statistical  nature  of  this  process. 

4  5  6 

The  papers  by  Plemons  and  Plemons,  Middleton,  and  Shooter  ’  contain 

experimental  results  of  the  temporal  and  spectral  characteristics  of  the 

process.  Their  experiment  utilized  a  single  transducer  to  measure  rever- 

7  8 

beration  signals  from  a  lake  surface.  In  a  similar  manner  Frazer  and 
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Omichi  studied  the  joint  spatial  and  temporal  properties  of  surface 
reverberation  by  analyzing  signals  from  multiple  receiving  elements. 
Although  these  works  have  contributed  significantly  to  our  knowledge  of 
this  particular  random  process,  there  is  still  a  large  gap  between  cur¬ 
rent  sonar  signal  processors  and  the  optimal  processors  mentioned  above. 

To  justify  the  implementation  of  more  powerful  processors,  the  sonar 
designer  needs  to  know  how  much  improvement'  there  will  be  over  simpler 
processors.  This  is  an  area  where  the  acoustic  signal  processing  litera¬ 
ture  is  not  very  abundant.  There  have  been  a  few  such  studies,  two  of 
which  are  mentioned  here.  Hof fman^ ’ ^  presented  some  theoretical  and 
experimental  work  implementing  a  simplified  (compared  to  Middleton's  and 
Groginsky's)  Bayesian  processor  using  sonar  echoes  from  several  contacts. 
He  used  the  processor  to  classify  echoes  from  objects  of  different  sizes 
and  materials.  To  learn  more  about  the  optimum  detector's  performance, 
there  is  a  need  for  experimental  evaluation  using  surface  reverberation 
data . 

12 

One  such  study  was  performed  by  Baker  and  Cunningham  using  single 
channel  beamformed  data  from  the  Navy's  AN/SQS-23  sonar.  These  aata 
included  submarine  echoes  as  well  as  segments  of  pure  surface  reverbera¬ 
tion.  Reverberation  and  signal  ensembles  were  formed  over  a  hundred  con¬ 
secutive  pings  and  used  to  compute  the  signal  and  noise  sample  covariance 
matrices.  Signal-plus-noise  returns  were  artificially  created  by  adding 
a  signal  return  to  the  reverberation  echoes  at  various  signal-to-noise 
ratios.  The  signals  thus  formed  were  processed  with  the  discrete  form 
of  the  Bayesian  detector  specified  by  Middleton  and  Groginskv. 


The  Present  Stud\ 


The  present  study  is  an  experimental  evaluation  of  an  optimum 
detector  for  a  multichannel  active  sonar  operating  in  the  presence  of 
surface  reverberation.  The  detector  is  optimum  in  the  Nevman-Pearson 
sense,  i.e.,  the  detection  prouability  is  maximized  lor  a  given  false 


alarm  probability.  The  design  of  the  detector  is  based  on  the  assumption 
of  a  known  signal  in  Gaussian  noise. 


Experiments  in  the  open  ocean  with  actual  Navy  sonars  are  expensive 
and  present  difficulties  in  controlling  many  parameters  of  interest.  On 
the  other  hand  a  high  degree  of  experimental  control  is  possible  in  taking 
lake  data,  and  lake  surface  reverberation  bears  many  similarities  to  that 
of  the  ocean.  In  this  study  an  ensemble  of  reverberation  returns  is 
created  by  scattering  sound  from  a  wind  driven  surface  of  a  freshwater 
lake.  Signals  on  each  of  six  receiver  elements  (5.1  cm  center  to  center) 
are  sampled  and  digitized.  Signal-plus-noise  returns  are  formed  in  a 
general  purpose  computer  by  adding  an  experimental  target  echo  to  each 
member  of  the  noise  ensemble  at  various  signal-to-noise  ratios.  The 
covariance  matrix  is  estimated  from  the  ensemble  and  employed  by  the 
detector.  The  detector's  performance  is  then  analyzed  by  processing  the 
signal-plus-noise  returns. 

The  detector's  performance  is  presented  in  the  form  of  receiver 
operating  characteristic  (ROC)  curves.  The  aspects  of  the  detector  in¬ 
cluded  in  this  study  are: 

1.  the  effect  of  sampling  frequency; 

2.  the  effect  of  transmit  type;  transmissions  include  the  following 

types:  a  1.0  msec  pui.sed  continuous  wave  (cw)  signal  at  a  carrier 

frequency  of  80  kHz  and  a  1.0  msec  pulsed  linear  frequency  modu¬ 
lated  (LFM)  signal  with  a  10  kHz  sweep  centered  at  80  kHz; 

3.  the  effect  of  assuming  the  noise  process  is  stationary  over  the 
time  duration  of  the  signal; 

4.  the  effect  of  assuming  the  noise  samples  are 

(1)  uncorrelated  in  space  but  correlated  in  time, 

(2)  uncorrelated  in  time  but  correlated  in  space, 

(3)  uncorrelated  in  both  space  and  time. 
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II.  THEORY 

This  section  includes  a  brief  theoretical  review  of  the  technique 
used  to  design  the  detector  being  considered.  The  Neyman-Pearson  lemma 
is  presented,  without  proof,  as  justification  for  choosing  the  likelihood 
ratio  test  employed  by  the  detector.  Several  assumptions  are  made  leading 
to  the  detector's  implementation  and  the  performance  prediction.  Finally, 
the  sampling  procedure  and  its  relation  to  the  detection  problem  are 
described . 

A.  Background 

In  an  abstract  sense  the  detection  problem  is  represented  by  Fig.  2.1. 
In  a  sonar  system  the  observation  might  be  a  set  of  samples  from  trans¬ 
ducer  element  signals  over  some  finite  time  interval.  The  observation 
will  be  considered  as  an  event  which  could  have  been  produced  by  one  of 
two  possible  causes.  The  two  possible  causes  considered  here  are  signal 
plus  noise  and  noise  only.  The  null  hypothesis,  that  noise  only  was 
responsible  for  the  event,  will  be  designated  as  Hq  and  the  other  hypoth¬ 
esis,  that  the  event  was  caused  by  signal-plus-noise,  by  H^. 

The  data  processor  must  take  an  observation  x  and  produce  a  decision 
as  to  which  hypothesis  is  most  likely.  The  decision  rule  will  be  deter¬ 
mined  by  dividing  the  observation  space  into  two  regions  Rq  and  so 
that  for  any  observation  x  either  Hq  or  can  be  chosen  as  the  true  case. 
To  arrive  at  a  suitable  decision  boundary  between  Rq  and  R^  the  four 
possible  outcomes  must  be  considered.  The  four  outcomes  are: 

(1)  given  is  true,  the  detector  decides  (true  alarm), 

(2)  given  Hq  is  true,  the  detector  decides  Hq  (true  dismissal) , 

(3)  given  is  true,  the  detector  decides  (false  alarm), 

(4)  given  is  true,  the  detector  decides  Hq  (false  dismissal). 
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If  the  probability  that  or  is  tine  is  known  a  priori  and  a 
cost  is  assigned  to  each  of  the  four  outcomes,  then  the  Bayes  criterion 
can  be  used  to  determine  a  decision  rule.  This  is  accomplished  by  mini¬ 
mizing  the  expected  cost  or  average  risk  R.  Let  nQ  and  be  the 
a  priori  probabilities  of  Hq  and  ,  respectively,  and  be  the  cost 
of  announcing  when  is  true.  Then  the  average  risk  R  is 

R  =  CLn  n  Pr{choosing  H  I H  is  true} 

00  o  o 1  o 

+  Cin  H  Pr{choosing  H  Ih  is  true} 

10  o  1 '  o 

+  nx  Pr{choosing  is  true} 

+  Pr{choosing  Hq  ( is  true} 


It  will  further  be  assumed  that  the  conditional  densities  of  the 
observation  x  under  both  hypotheses  are  known.  Letting  p^(x)  and  p^(x) 
represent  the  conditional  densities  of  x  under  the  null  and  true  hypoth¬ 
eses,  respectively,  the  risk  is 


/ 


R  =  C00  no  /  PoU)dx  +  C10  Fo  f  Po(x)dx 


+  C11  n 


1  J  p]/x)dx  +  coi  ni  /  Pj (x) dx 


With  some  manipulation  it  can  be  shown  that  choosing  the  decision 
boundary  in  this  manner  results  in  including  those  values  of  x  for 
which 


P1(x)  (C 

Pq(x)  1^(0 


TO 

01 


-c 


00 

11 


) 

) 


The  quantity  on  the  left  is  called  the  likelihood  ratio  and  is  designated 
by  A(x).  The  quantity  on  the  right  represents  a  threshold  and  is  labeled 
A^,.  This  test  is  called  a  likelihood  ratio  test  and  can  be  simply  written 
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A(x)  <  A 
H 

o 

The  implementation  of  this  test  is  shown  in  Fig.  2.2. 

B.  The  Neyman-Pearson  Lemma 

In  many  cases  neither  the  a  priori  probabilities  and  IIo  nor  the 
costs  of  C.,  are  known.  In  this  study  including  these  costs  and  prob¬ 
abilities  would  only  complicate  the  evaluation  of  the  optimal  data  pro¬ 
cessor.  Therefore,  a  more  general  criterion  is  needed  for  the  formulation 
of  the  detector  considered  here.  The  criterion  is  provided  by  Neyraan 
and  Pearson  and  is  based  on  a  theorem  discussed  below. 

Before  proceeding  with  the  so-called  Neyman-Pearson  lemma,  a  few 
statistical  terms  related  to  hypothesis  testing  are  introduced.  The 
first  of  these  terms  is  an  error  of  the  first  kind.  This  refers  to  the 
third  outcome  of  the  decision  process  mentioned  above  (saying  when 
is  true).  The  conditional  probability  associated  with  this  error  is 
called  the  false  alarm  probability  or  the  size  or  level  of  the  test  and 
is  conventionally  designated  by  a.  The  conditional  probability  of  a 
true  alarm  (saying  when  is  true)  is  called  the  probability  of  detec¬ 
tion,  or  the  power  of  the  test,  and  is  generally  designated  by  6.  The 
Neyman-Pearson  criterion  consists  of  maximizing  the  power  of  the  test  6 
with  the  size  of  the  test  a  constrained  to  be  less  than  or  equal  to  some 
prescribed  value. 

The  following  presentation  of  the  Neyman-Pearson  lemma  is  taken  from 
Lehmann. ^ 

Let  <J)(x)  be  the  probability  of  announcing  when  the  observation  is 
x.  In  the  context  of  the  Bayes  test,  4>(x)  would  be  either  zero  or  one, 
depending  on  what  region  x  lies  within.  Let  E^f^tx)}  be  the  expected 
value  of  <j)(x)  under  the  null  hypothesis.  In  addition,  let  P  and  be 


10 


ANNOUNCE  H, 


FIGURE  2.2 

LIKELIHOOD  RATIO  PROCESSOR 


the  probability  distributions  possessing  derivatives  pQ  and  p^, 
respectively,  with  a  measure  p. 

(i)  Existence.  For  testing  HQ:po  against  the  alternative  H^:p^, 
there  exists  a  test  <p  and  a  constant  A such  that 

Eq { 4>  (x)  }  =  a  (2.1) 

and 

1  when  p1(x)  >  AtPq(x) 

0  when  p^(x)  <  A^p^(x) 

(ii)  Sufficient  conditions  for  a  most  powerful  test.  If  a  test 

satisfies  Eq.  (2.1)  and  Eq .  (2.2)  for  some  A^,  then  it  is  most  powerful 

(maximum  B)  for  testing  pQ  against  at  level  a. 

(iii)  Necessary  condition  for  a  most  powerful  test.  If  4>  is  most 

powerful  at  level  a  for  testing  pQ  against  p^,  then  for  some  it  sat¬ 

isfies  Eq.  (2.2)  almost  everywhere  with  respect  to  the  measure  p.  It  also 
satisfies  Eq.  (2.1)  unless  there  exists  a  test  of  size  less  than  a  and  with 
power  equal  to  one. 


(2.2) 


The  test  4>(x)  in  part  (i)  means  that,  given  the  observation  x, 
hypothesis  is  announced  with  probability  one  if  (p^ (x) /pQ (x) )  is 
greater  than  A  .  Therefore,  <J> ( x)  is  a  likelihood  ratio  test  with  thres- 
hol  A  .  The  threshold  here  is  not  a  function  of  costs  or  a  priori 
probabii ities  as  in  the  Bayes  test,  but  rather  it  is  a  function  of  the 
conditional  densities  pQ  and  p  and  the  size  of  the  test  a.  To  see  this, 
consider  Eq .  (2.1)  where 
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A^  is  chosen  so  that  the  value  of  the  integral  just  equals  a.  Note  that, 
for  At=0,  the  integral  is  over  all  x  and  is  equal  to  unity.  The  power 
of  the  test  at  level  a  is  given  by  the  expected  value  of  <J>(x)  under  the 
hypothesis  as  follows. 


8  =  E1{*(x)> 


4>(x)p1(x)dx  = 


P1(x)dx 


All  x 


'x: 


P1(x ) 

PG(x)  > 


(2.4) 


Summarizing  the  test,  it  can  be  said  that,  given  a  false  alarm  probability 
a,  a  threshold  for  the  likelihood  ratio  is  chosen  by  Eq .  (2.3).  Given 
this  A^  the  detection  probability  8  is  computed  from  Eq.  (2.4). 

An  evaluation  of  a  given  processor  (i.e.,  given  p^  and  p^)  can  then 
be  made  by  considering  a  plot  of  8  versus  a.  This  plot  is  the  so-called 
receiver  operating  characteristic  (ROC)  curve  and  is  utilized  herein  to 
present  the  detector's  performance.  The  thrust  of  the  Neyman-Pearson 
lemma  is  that,  given  p^  and  p-^,  there  exists  no  test  whose  ROC  curve  is 
above  the  curve  associated  with  4> (x) ,  the  likelihood  ratio  test.  If 
there  were  such  a  test  it  would  be  more  powerful  at  the  level  a  where 
the  curve  was  greater  and  hence  would  be  a  contradiction  of  the  lemma. 

C.  The  Likelihood  Ratio  Detector 


This  section  gives  the  assumptions  upon  which  the  detector  is  based. 
The  specific  detector  structure  is  formulated  and  the  method  for  perform¬ 
ance  evaluation  is  given. 

1 .  Assumptions 

The  major  assumptions  leading  to  the  detector  being  studied 
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here  are  as  follows: 


(1)  Detection  is  a  binary  operation,  i.e.,  one  of  two  hypotheses 
must  be  chosen  given  an  observation:  Hq:  observation  due  to  noise  alone, 
and  H^:  observation  due  to  signal-plus-noise.  Only  one  observation  is 
made  even  though  in  many  practical  sonar  problems  multiple  "looks"  at  a 
potential  target  are  possible  before  the  detection  decision  is  reached. 

(2)  Processing  is  done  on  a  finite  number  of  samples  taken  during 
a  fixed  time  interval.  This  interval  corresponds  to  the  length  of  the 
known  signal  which  is  approximately  equal  to  the  pulse  length  of  the 
transmission.  In  the  sonar  detection  problem,  the  optimal  detector  would 
be  time  v  irving  since  the  surface  reverberation  process  is  considered  to 

g 

be  nonstationarv.  In  other  words  an  actual  system  would  process  several 
fixed  time  intervals  each  corresponding  to  a  different  point  in  range. 

(3)  The  signal  is  known,  consisting  of  samples  taken  from  an  echo 
from  a  submerged  styrofoam  ball.  Therefore,  range  and  bearing  to  the  tar¬ 
get  (when  present)  are  known  as  well  as  the  shape  of  the  echo.  The  more 

general  problem  of  a  random  signal  with  unknown  target  location  is  beyond 

3 

the  scope  of  this  paper.  As  Middleton  and  Groginsky  point  out,  the 
assumption  of  known  location  is  not  purely  academic.  A  series  of  such 
optimal  detections  might  be  carried  out  under  this  constraint  at  differ¬ 
ent  ranges  and  bearings. 

(4)  Both  system  and  ambient  noise  are  considered  to  be  negligible 
in  relation  to  the  noise  level  generated  by  surface  reverberation.  This 
is  generally  the  case  in  most  active  sonar  problems  at  close  ranges. 

(5)  The  noise  process  is  Gaussian  but  not  necessarily  zero  mean  or 
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stationary.  The  work  by  Frazer  ’  indicated  that  this  assumption  is  valid 
on  a  first  order  basis,  i.e.,  the  first  order  distribution  function  is 
Gaussian.  The  assumption  that  the  process  is  jointly  normal  both  in  time 
and  space  is  not  quite  as  well  substantiated.  However,  if  only  informa¬ 
tion  about  the  second  order  moments  is  available,  one  may  reasonably 
assume  the  process  is  Gaussian.  This  statement  is  based  on  the  fact  that, 
given  an  autocorrelation  function  of  a  random  process,  there  always  exists 
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a  Gaussian  random  process  possessing  the  same  autocorrelation  function. 
Therefore,  the  actual  and  the  assumed  processes  are  indistinguishable 
from  the  standpoint  of  what  is  known. 


(6)  The  covariance  of  the  noise  process  is  known.  In  general  this 

is  not  easy  to  determine;  however,  it  is  possible  that  the  covariance 
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function  might  be  predicted.  For  the  optimum  detector  considered  here, 

the  covariance  function  is  assumed  known  precisely. 

(7)  Signal  and  noise  are  additive.  Since  the  noise  considered  in 
this  study  results  from  reverberation  of  sound  waves  from  a  rough  surface, 
it  is  reasonable  that  these  sound  waves  would  add  to  an  echo  from  a  tar¬ 
get  at  approximately  the  same  range  from  the  receiver  as  the  surface, 

2.  Detector  Formulation 


In  the  previous  development  of  the  likelihood  ratio  detector, 
the  observation  x  could  have  been  a  vector  as  well  as  a  scalar  quantity. 
Such  is  the  case  for  the  input  to  the  data  processor  considered  here. 

The  development  of  this  processor  proceeds  as  follows.  Consider  M  sets 
of  received  data. 


V(t)  =  [v-1(t),i’2(t),...,VM(t)]T 
T 

where  []  indicates  transpose  operation,  and  V^(t)  is  the  waveform  observed 
at  the  ith  element  over  a  fixed  time  interval.  They  are  assumed  to  have 
been  generated  by  either  ikCt)  or  6i(t)+ui(t)  corresponding,  respectively, 
to  the  hypothesis  Hq  that  noise  alone  is  present  or  to  ,  that  signal- 
plus-noise  is  present.  Next  the  v\(t)'s  are  sampled  in  time,  resulting 
in  a  set  of  M  column  vectors, 


V 


Iv1,v2,...,VmJ 


where  each  _v^ 


is  the  column  vector. 


vi  =  ’ '  •  •  ,yi/tN^T’  1  =  1,M 

and  N  is  the  number  of  samples  taken  during  the  fixed  time  interval.  In 
the  sampled  data  problem  the  two  hypotheses  are  then 

H  V  -  N 

o  —  — 


and 


H  :  l  =  S  +  w 


where 


isi  [ n^  ,  1^2  *  •  -  1 3  ’ 

T 

-i  =  *-Ui ( t1)  ’,Zi(t2) . Hi(tN^  ’  i  =  l.M 

S  =  [-l’-2’  •  •  •  ’ 


and 


-i  =  f5i(t1)>ii(t2)’*--’5i(tN) 1  ’  i  =  1,M 


First  the  likelihood  ratio  is  formed  by 


A(i/)  =  - 


PG(^) 


(2.5) 


where  p  ((/)  and  pQ(j/)  are  the  joint  density  functions  of  under  and 

H  ,  respectively.  Under  H  ,  l'=N  and  N  is  assumed  Gaussian;  therefore, 
o  o -  — 


Po(-}  =  ^ 


j-  l/2(2-^)T  k'V-M)J 
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where 


M  =  noise  mean  vector 
=  E{W}, 

K  =  NM  x  NM  positive-definite,  symmetric  covariance  matrix 

=  e{  (W-M)  (w-M)t}.  and 


Under  H  ,  and  the  mean  vector  and  covariance  matrix  are 

E{U}  =  E{S+M)  =  s  +  M 


and 


E^(U-E{1/})  (V-E{V})T}  =  Ej[S  +  M  -  (S+M)  ]  \S  +  M  -  (S+M)]T( 

=  E  j(N-M)  (M-M ) T  | 

=  K 

Therefore, 

p^(l/)  =  exp{-l/2(l/-M-S)^  tC  (U-M-5) } 

Inserting  these  expressions  for  p  and  p^  into  Eq .  (2.5),  simplifying, 
and  taking  the  logarithm  yields 

InA(U)  =  1/Tk_  15  -  MTK_1S  -  1/2(SV1S) 

The  likelihood  ratio  test  is 

H1 

lnA(U)  In  A? 

H 
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or,  equivalently. 


X  ^  j/T  K_1  S  <  InA  +  MT  K-1  S  +  1/2 (STK  1S)  =  X 

H  T  1 

o 

This  last  step  is  possible  because  the  term  on  the  left  is  the  only  one 
which  depends  on  the  observation  (/.  The  other  terms  are  functions  of 
known  quantities  and  consequently  can  be  considered  as  part  of  the  thres¬ 
hold  which  will  be  determined  later. 

Before  proceeding  further  with  the  analysis,  a  few  words  are  in 
order  concerning  the  nature  of  the  above  data  processor.  Letting 

6  =  K_1S  ,  (2.6) 

the  output  X  can  be  written  as  the  scalar  product  of  and  B^, 

X  =  l/TB  .  (2.7) 


Also  let 


where  the  b.'s  represent  the  filter  coefficients  for  channel  i.  A  block 


diagram  of  the  multichannel  likelihood  ratio  detector  is  shown  in  Fig.  2.3 


Consider  K  as  a  matrix  of  submatrices  C..  each  of  order  N  (the 


number  of  time  samples)  as  follows: 


,-l 


kl 

—12  ’ ’ ’ 

-1M 

c 

c  „  ... 

r 

-21 

—22 

-2M 

c 

c 

c’ 

-Ml 
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FIGURE  2.3 

MULTICHANNEL  LIKELIHOOD  RATIO  PROCESSOR 
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This  notation  is  not  meant  to  imply  that  C. .  is  only  a  function  of  the 
covariance  matrix  between  channels  i  and  j.  On  the  contrary,  since  C_ 
represents  a  subsection  of  the  inverse  of  the  complete  covariance  matrix 
I<,  each  of  the  elements  of  (?  is  a  function  of  al  1  the  elements  of  K. 
Using  this  notation  the  elements  h.^  can  be  expressed  as 


M 

h =  E  L:ik^’  1  =  1*M  • 

k  =  l 

It  is  interesting  to  note  that  in  general  the  b^'s  depend  not  only  on  the 
signal  and  noise  at  channel  i,  but  also  on  the  signals  and  noise  on  each 
of  the  other  channels  as  well. 


Further  insight  into  the  detector  structure  may  be  gained  by 
considering  a  few  special  cases.  First  consider  a  stationary  Gaussian 
process  with  unit  variance  which  is  independent  both  in  time  and  space. 
The  inverse  covariance  matrix  in  that  case  is  the  identity  matrix  and  the 
detector  is  a  set  of  filters  _b . ,  i=l,M,  each  matched  to  the  expected  sig¬ 
nal  at  channel  i.  The  detector  output  1  in  that  case  is 


\ 


i  =  l 


‘•ftwv  V  • 

i=l  j=l 

The  role  of  the  covariance  matrix  may  be  observed  by  relaxing 

the  assumption  of  stationarity .  Thus,  if  the  noise  variance  changes  with 

time  and  channel,  then  the  diagonal  elements  of  the  inverse  covariance 

->  2 

matrix  are  given  by  l/o..,  where  a . .  is  the  variance  ot  the  j th  time  sam- 

lj  tj 

pie  on  channel  i.  Consequently  the  vector  of  filter  coefficients  for 
channel  i  is 
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b. 
— 1 


W 

2 

il 


w 

2 


ii(tN) 


a  -vt 

lN 


It  Is  apparent  that  the  test  statistic  consists  of  weighting  the  input 
.  ? 

samples  v. (t.)  by  5.(t.)/a..  as  follows 
i  j  7  i  J  ij 


X 


EE 

i=l  j  =  l 


(cj 


Vhi 

o 

o7 . 

ij 


In  other  words,  if  the 

expected  to  be  large  in 

then  the  sample  i'.(t.) 

1  J 


signal  is  present  and  a 
absolute  value  compared 
is  weighted  accordingly. 


particular  sample  is 
to  the  noise  variance 


Another  interesting  case  arises  when,  in  addition  to  the 
previous  assumptions  of  stationarity ,  independence,  and  normality,  it  is 
also  assumed  that  the  signal  is  identical  on  each  channel.  (This  condi¬ 
tion  could  be  created  by  adding  a  suitable  time  delay  and  weight  factor 
in  each  channel,  assuming  the  signal  shape  is  identical  on  all  channels.) 
Then  the  test  statistic  becomes 


A 


and  the  processor  amounts  to  a  summation  of  the  inputs  at  each  channel 
followed  by  a  single  matched  filter.  This  form  of  signal  processing  is 
prevalent  in  modern  sonars. 

3.  Threshold  Determination  and  Performance  Prediction 


In  accordance  with  the  Noyman-Pearson  criterion,  the  threshold 
for  the  likelihood  ratio  test  is  chosen  to  achieve  a  given  false  alarm 
probability,  a.  In  general  an  analytic  expression  for  the  threshold 
as  a  function  of  a  is  unobtainable  and  the  relation  must  be  determined 
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experimentally.  Since  the  noise  is  assumed  Gaussian  in  this  study,  the 
functional  relationship  can  easily  be  derived. 

r  _  i 

If  \-L  K  S',  then  it  represents  a  linear  combination  of  the 
elements  of  \J_,  assumed  Gaussian;  therefore,  A  is  itself  a  normal  random 
variable.  Its  expected  value  can  be  written 


K { X }  =  E{VTK-1S}  =  E{l/T}K 


The  variance  of  \  it 


Va  r f ' ^  = 


=  =  E | [ A-E{ A } ] 2| 

=  k{[  (^-E{I/1)TK_1S)]2|- 

=  E-jsT(K_1)T(^-E{^})  (V-E{l'))V1sj 

=  5T(K_1)T  E{(V-E{j/})  0/-E{V})TK-1sJ. 


=  St(k  V  K  K_1S 


„ .  ,  .-1  .  T  -1 

Since  n  is  symnictru:,  (K  )  =K  and 


2  CT  -lc 
o  ^  =  6  h  o 


(2. 81 


The  threshold  and  the  resulting  power  of  the  test  are  determined  by  con¬ 
sidering  the  distribution  of  A  under  the  null  and  true  hypotheses.  Under 

hvpotbos  i  s  II  ,  A’s  mean  value  is 
o 


!■:  {Al  =  m  =  E  {I/T}K_15 
o  o  o  —  —  — 


III 

o 


5 


>  ■> 


or 


Under  hypothesis  H^,  the  mean  of  X  is 


EL{  A  >  =  irij  =  E1{UT}K~1S 

=  (MT+ST)K-1S 


or 


T  -1  T  -1  ? 

m  =  M  K  S  +  S  K  S  =  m  +  oT 

1 - —  —  -  —  o  \ 


(2.9) 


As  was  demonstrated  earlier,  the  variance  of  X  is  the  same  under  both 
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hypotheses  and  is  represented  hy  a~.  Thus  it  is  apparent  that  the  data 
processor  reduces  the  samples  of  the  vector  random  process  1/  to  a  sample 
of  a  normal  random  variable  X.  Typical  conditional  probability  density 
functions  for  X  given  Hq  and  are  shown  in  Fig.  2.4  as 
P,/H(x).  respectively. 

Since  the  likelihood  ratio  test  can  be  represented  by 


P^/._(A)  and 
o 


X  <  X 

H 

o 

the  false  alarm  probability  is  given  by 

CIO  CO 

«.  -f  P1/H(*>«  -j  ^  <*r 
X.  X 

1  I 

This  probability  is  indicated  by  the  shaded  area  in  Fig.  2.4(a). 
Making  the  change  of  variables. 
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(a)  Ot  =  FALSE  ALARM  PROBABILITY 


(b)  0  =  DETECTION  PROBABILITY 
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FIGURE  2.4 

TYPICAL  CONDITIONAL  PROBABILITY  DENSITIES  FOR  THE  TEST 
STATISTIC  X  AND  ASSOCIATED  ERROR  PROBABILITIES 


(? .  in 


and  finally. 


XT  =  mQ  +  <J)-1  (1-a)  .  (2.10) 

Note  that  if  the  noise  process  is  zero  mean,  then  M=()  and  consequently 

m  =0. 

o 


Given  the  threshold  X^,  the  resulting  detection  probability  0 
can  be  predicted  by  considering  the  density  function  for  X  under  the 
true  hypothesis.  Thus  the  expression  for  0  is 


6 


This  probability  is  indicated  by  the  shaded  area  in  Fig.  2.4(b).  By 
the  appropriate  change  of  variables, 


a  f  1  -U2/2  , 

'■j  S'  du 


XT_m1 


XT~ml 


A 

-  f  -L  .V 
J  Jli 


/  2 


du 


6  =  1-$ 


X,  -  mx' 


(2.11) 


An  expression  for  0  as  a  function  of  a  may  be  obtained  by 
substituting  Eqs.  (2.9)  and  (2.10)  in  F.q.  (2.11).  Therefore, 
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e  =  1  -  $ 


=  1-0 


r  -i .  "i 

K +  aA $  (1_a)J 

-  mx 

a. 

A 

m  +  a.  <i  ^  ( 1-a)  - 

m  - 

2l 

o  A 

0 

X 

or 


6  =  1  -  4>|o  1  ( 1  -a )  -  a^j  .  (2.12) 

This  is  Che  predicted  ROC  curve  for  the  detector  and  is 

presented  in  Fig.  2.5  for  several  values  of  o.  .  Horton‘S  had  demonstrated 

that  a ^  is  related  to  the  output  signal-to-noise  ratio.  Intuitively  this 

can  be  seen  by  observing  that  the  separation  between  the  two  curves  in 
2 

Fig.  2.4  is  o^.  Since  the  detector's  performance  improves  with  increasing 

a^,  one  might  be  led  to  consider  Eq.  (2.12)  in  that  regard.  If  the  noise 

were  independent  of  the  signal  then  S  could  be  chosen  to  maximize  o  . 

A 

However,  in  the  active  sonar  problem  the  noise  is  very  definitely  related 
to  S  since  both  are  produced  by  the  transmitted  pulse.  Thus,  any  change 
in  -S  caused  by  altering  the  transmission  results  in  a  different  covariance 
matrix. 


One  of  the  objectives  of  this  study  is  to  compare  the  predicted 
ROC  curves  with  experimental  results  obtained  by  processing  experimentally 
acquired  acoustic  data  with  the  proposed  detector.  The  manner  in  which 
the  experimental  ROC  curves  are  derived  will  be  described  in  Section  IV. 

D.  Sampling  Technique 

It  was  stated  earlier  that  the  input  to  the  data  processor  would  be 
a  set  of  samples  from  a  fixed  time  interval.  Since  the  processor's  compu¬ 
tational  complexity  is  roughly  proportional  to  the  square  of  the  number  of 
samples  (due  to  the  matrix  inversion  of  K) ,  the  sampling  frequency  and 
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techniques  are  major  concerns.  In  this  section  the  sampling  method  is 
described. 


In  many  sonar  applications  the  signals  of  interest  are  bandlimited 
to  the  extent  that  most  of  the  energy  is  concentrated  within  a  band  of 
frequencies.  This  is  generally  the  case  with  an  active  sonar  which  trans 
mits  acoustic  energy  at  a  bandwidth  which  is  small  compared  to  the  car¬ 
rier  frequency.  Also,  the  first  stage  of  signal  processing  by  the 
receiver  usually  involves  bandpass  filtering.  Such  is  the  case  for  the 
signals  considered  here  for  processing  by  the  likelihood  ratio  detector. 

To  facilitate  the  discussion  of  the  sampling  technique,  the  signal 
is  assumed  to  be  a  real  function  of  time  which  is  strictly  bandlimited  to 
the  frequency  interval  (fQ  -  W/2)  <  [f[  <  (fQ  +  W/2).  Furthermore,  it  is 
assumed  that  the  bandwidth  W  is  small  compared  to  the  carrier  frequency 
fQ.  Such  a  signal  can  be  represented  in  the  time  domain  by 

v(t)  =  A(t)  cos [ 2rr f ^ t-+-<J>  ( t )  ]  ,  (2.13) 

where  A(t)  and  <}>(t)  are  the  amplitude  and  phase  functions,  respectively. 

The  quadrature  representation  obtained  by  the  identity 
cos(A+B)  =  rosA  cosB  -  sinA  sinB  is 

v’(t)  =  X(t)  cos(2iTfot)  +  Y(t)  sin(27rfot)  , 


where 


X(t)  =  A(t)  cosiJ>(t),  and 
Y(t)  =  -A(t)  sin<f>(t) 


X(t)  and  Y(t)  are  called,  respectively,  the  in  phase  and  out  of  phase 
quadrature  components  of  v/(t)  and  are  bandlimited  to  the  interval  |f| 


<W/2 
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X(t)  and  Y(t)  can  be  completely  represented  by  their  samples  taken  at  a 

minimum  rate  of  W  samples  per  second,  commonly  referred  to  as  the  Nyquist 
16 

rate.  » 

One  method  of  obtaining  samples  of  X(t)  and  Y(t)  consists  of 

heterodyning  u(t)  appropriately  to  produce  the  two  quadrature  components 

and  then  sampling.  A  preferred  method  is  to  obtain  samples  of  X(t)  and 

Y(t)  by  sampling  u ( t )  directly  in  the  manner  presented  by  Pitt  and 
17-19 

Grace.  A  simplified  version  of  the  technique  is  described  here. 

Samples  of  X(t)  are  obtained  by  sampling  v(t)  at  the  times 
t=nk/fQ,  n=l,2,3, . . . ,N/2,  where  the  sampling  parameter  k  is  a  positive 
integer  which  determines  the  time  interval  between  samples  of  X(t)  in 
multiples  of  1 / f ^ ,  and  N  is  the  total  number  of  samples  from  both  com¬ 
ponents.  To  ensure  that  samples  of  X(t)  are  taken  at  a  rate  no  less  than 
the  bandwidth  W,  k  must  be  less  than  or  equal  to  the  largest  integer  that 
is  less  than  or  equal  to  fo/W.  If  v(t)  is  also  sampled  at  the  same  fre¬ 
quency  but  at  times  t= (nk/fQ)+(l/4fo) ,  n=l , 2 , 3 , . . . ,N/2 ,  samples  of  Y(t) 
are  produced.  These  relationships  are  obtained  by  substituting  the  above 
times  in  Eq.  (2.13).  For  more  details  the  reader  is  referred  to 
Refs.  17-19. 


Oftentimes  the  amplitude  and  phase  functions  of  v(t)  are  of  interest. 
They  may  be  computed  from  samples  of  X(t)  and  Y(t)  as  follows: 

A(t)  =  VX2(t)  +  Y2  (t) 


and 


<J>  (c) 


tan 


-1 


-Y  ( t) 
X(t) 


However,  since  the  samples  of  Y(t)  occur  a  quarter  cycle  of  the  carrier 

later  (l/4f  )  than  those  of  X(t) ,  the  following  approximation  must  be 
o 

made . 
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If  W  is  small  compared  to  f^,  then  Y(t)  is  slowly  varying  with  respect 
to  cos(2irfot)  and  the  approximation  is  a  good  one. 

To  study  the  effect  of  sampling  frequency  on  the  detector's 
performance  the  quadrature  sampling  technique  is  employed.  An  input 
vector  v^  from  a  given  channel  i  consists  of  quadrature  samples  of  v . (t) 
as  follows: 


-i  =  (ui(t:i>>  Ui(tl>’  ui(t2)>  Ui(tN/2)’  Ui(tN/2)]  ’  (2,14) 


where 


t 


n 


corresponding  to  the  in  phase  component,  and 


t' 

n 


(2.15) 


(2.16) 


yielding  samples  of  the  out  of  phase  component.  The  effect  of  sampling 
frequency  on  detector  performance  is  observed  by  varying  the  sampling 
parameter  k. 


III.  EXPERIMENTAL  PROCEDURE 


This  chapter  describes  the  acquisition  of  data  used  in  the  evaluation 
of  the  likelihood  ratio  detector.  These  data  consist  of  measurements  of 
reverberation  from  the  lake  surface  and  echoes  from  a  styrofoam  sphere. 

The  reverberation  data  were  obtained  in  an  experiment  conducted  in 
September  1974.  This  experiment  is  described  briefly  here  and  in  greater 
detail  in  Ref.  9.  A  second  experiment  was  conducted  in  July  1975,  in 
which  the  styrofoam  sphere  echoes  were  recorded  using  the  same  equipment 
as  the  1974  experiment.  A  notable  improvement  relating  to  the  alignment 
of  data  from  different  channels  (see  Section  IV. A)  was  made  in  the  sec¬ 
ond  experiment.  The  second  data  collection  exercise  was  also  conducted 
to  gather  additional  reverberation  data;  however,  these  data  were  deemed 
unsuitable  for  this  study  due  to  a  low  interchannel  coherence.  A 
description  of  the  arrangement  and  location  of  equipment  is  given 
including  types  of  sonar  transmissions  utilized.  Finally,  the  analog 
recording  and  analog-to-digital  conversion  processes  are  discussed. 

A.  Data  Acquisition  Site  and  Environmental  Conditions 

The  data  used  in  this  study  were  obtained  at  Applied  Research 

20 

Laboratories  Lake  Travis  Test  Station  (LTTS) .  Figure  3.1,  adapted  from 
Ref.  20,  indicates  the  geographic  location  of  LTTS.  Since  LTTS  is  on 
the  south  end  of  the  lake,  data  are  taken  with  the  wind  from  a  northerly 
direction.  This  produces  a  fairly  rough  surface  so  that  meaningful 
reverberation  phenomena  may  be  observed.  Table  I  presents  the  pertinent 
environmental  factors  observed  during  the  reverberation  and  target  echo 
experiments . 
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FIGURE  3.1 

LOCATION  OF  LAKE  TRAVIS  TEST  STATION 


ENVIRONMENTAL  FACTORS  OF  REVERBERATION  AND 
TARGET  ECHO  EXPERIMENTS 


The  transducer  array  consisted  of  a  linear  receiving  array  with  a 
projector  rigidly  attached  (see  Fig.  3.2,  adapted  from  Ref.  9).  The 
receiving  array  was  constructed  of  20  rectangular  elements  with  planar 
surfaces.  Each  element  was  10.2  cm  high  and  4.95  cm  wide.  The  horizon¬ 
tal  separation  (edge  to  edge)  between  elements  was  0.15  cm.  Only  the 
center  six  elements  were  used  in  these  experiments.  The  projector  was 
also  planar  and  rectangular  with  a  horizontal  dimension  of  6.6  cm  and  a 
vertical  dimension  of  38.1  cm. 

Beam  patterns  (intensity  versus  angle)  for  receiver  elements  and 
projector  are  documented  in  Ref.  9.  The  beam  patterns  were  not  measured 
during  the  second  experiment  (July  1975)  ;  however,  they  are  assumed  to 
be  comparable  to  the  patterns  measured  in  the  1974  experiment. 

C.  Reverberation  Measurement 

The  geometry  for  the  reverberation  experiment  is  shown  in  Fig.  3.3 
(adapted  from  Ref.  9) .  A  transducer  array  was  placed  at  a  depth  of  3.6  m 
and  tilted  approximately  16°  toward  the  surface.  Two  transmit  signals 
were  employed: 

(1)  1.0  msec  pulsed  80  kHz  continuous  wave  (cw)  signal,  and 

(2)  1.0  msec  pulsed  Linear  frequency  modulated  (LF>1)  signal  with  a 
10  kHz  sweep  centered  at  80  kHz. 

A  total  of  99  reverberation  returns  were  generated  for  the  cw  data  and 
103  for  the  LFM  data.  The  time  separation  between  consecutive  transmis¬ 
sions  was  0.6  sec. 


Target  Echoes 


The  target  utilized  in  this  study  consisted  of  several  styrofoam 
disks  held  together  by  metal  plates  (see  Fig.  3.4).  The  overall  dimen¬ 
sions  are  approximately  30.5  cm  diameter  and  30.5  cm  height.  Echoes  from 
this  reflector  were  obtained  by  tilting  the  transducer  array  below  the 
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A  FRONT  VIEW  OF  THE  TRANSDUCER  ARRAY 
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horizontal  so  that  the  projector  axis  pointed  toward  the  target,  which 
was  located  at  a  horizontal  range  of  60  m  and  submerged  to  a  depth  of  9  m. 
This  depression  angle  was  chosen  to  minimize  the  surface  reverberation 
level  so  that  high  signal-to-background  echoes  might  be  obtained.  The 
transmit  signals  employed  here  were  the  same  as  those  used  in  the  rever¬ 
beration  experiment. 

E.  Data  Acquisition 

In  this  section  the  data  acquisition  system  employed  in  the  July 
1975  experiment  is  described.  The  earlier  experiment  was  similar  and 
is  described  in  Ref.  9. 

A  block  diagram  of  the  recording  system  is  shown  in  Fig.  3.5.  As 
shown  in  the  diagram,  several  signals  are  summed  with  the  received  sig¬ 
nal  before  being  led  to  the  tape  recorder.  These  signals  are  included 
to  facilitate  the  analog-to-digital  (A/D)  conversion  process.  A  block 
diagram  of  the  A/D  system  is  shown  in  Fig.  3.6,  adapted  from  Ref.  9. 

A  12-bit  A/D  is  used  coverin’  a  range  of  ±10  V. 

A  clock  signal  was  added  to  each  channel  to  serve  as  a  sample  clock 
for  the  A/D  converter.  Providing  this  signal  on  every  channel  permits 
later  compensation  for  the  pf<rects  f  t  the  tape  recorder's  dynamic  skew 
on  the  relative  phase  relationships  between  channels.  The  A/D  samples 
the  data  channel  at  every  zero  crossing  of  the  160  kHz  clock  signal, 
thus  yielding  a  uniform  sampling  frequency  of  320  kHz  (9fQ)  . 

The  sync  signal  was  added  to  provide  a  start  pulse  for  the  A/D 
converter.  As  seen  in  Fig.  3.6,  the  A/D  was  started  by  triggering  an 
oscilloscope  which  in  turn  provided  an  A/D  start  signal  after  an  adjust¬ 
able  delay. 

Since  each  tape  recorder  channel  was  played  back  and  digitized 
separately,  a  method  is  required  for  ensuring  that  the  time  base  is  the 
same  for  all  channels.  This  is  accomplished  by  recording  a  Barker  coded 

AO 


J _ J  I _ L. _ I 


FIGURE  3.5 

A  BLOCK  DIAGRAM  OF  THE  ANALOG  DATA  RECORDING  SYSTEM 
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FIGURE  3.6 

BLOCK  DIAGRAM  OF  DIGITIZATION  EQUIPMENT 
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pulse  at  80  kHz  center  frequency  shortly  after  the  transmit  signal. 

This  pulse  consisted  of  a  16-bit  binary  code  with  each  "1"  bit  repre¬ 
sented  by  two  cycles  of  the  80  kHz  carrier  at  one  phase  angle  and  each 
"0"  bit  by  two  cycles  180°  out  of  phase  with  the  "1"  bits.  The  particu¬ 
lar  code  employed  was  0001010110011111.  By  correlating  the  Barker  pulse 
on  a  reference  channel  and  ping  with  that  from  every  other  channel  and 
ping,  the  error  in  the  time  base  is  minimized  (see  Section  IV. A). 

Time  base  error  is  further  reduced  by  also  recording  a  reference  cw 
pulse  at  80  kHz  immediately  following  the  Barker  pulse.  If  the  time  base 
is  correct  for  all  pings  and  channels,  then  the  constant  phase  of  the 
reference  pulse  will  be  the  same  in  all  cases.  Using  this  reference 
pulse,  data  may  be  corrected  for  any  difference  found.  While  the  above 
procedure  is  not  foolproof,  the  results  (described  in  Section  IV. A)  indi¬ 
cate  the  method  is  highly  satisfactory. 

The  data  collection  procedure  was  very  similar  in  the  1974 
experiment  except  no  Barker  pulse  was  recorded. 
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IV.  EXPERIMENTAL  RESULTS 


The  results  of  the  experimental  performance  analysis  are  presented 
in  this  section.  First,  the  computational  procedure  for  aligning  the 
data  in  time  is  discussed.  The  computation  of  the  covariance  matrix 
estimate  used  by  the  detector  is  included  as  well  as  plots  of  the  cw  and 
LFM  matrices.  Details  of  the  detector  evaluation  are  presented  including 
the  computation  of  experimental  ROC  curves.  The  analysis  of  results 
considers  the  detector  performance  for  various  parameters  of  interest 
and  concludes  with  a  comparison  of  the  optimum  detector  to  several 
simplified  processors. 

A.  Alignment  of  Digital  Samples 

The  first  area  of  experimental  results  to  be  discussed  concerns  the 
alignment  of  data  from  different  pings  and  different  channels  of  the 
recording  system.  The  method  consists  of  time  shifting  data  from  one 
event  with  respect  to  another  by  an  amount  which  yields  the  greatest 
crosscorrelation  of  the  two  Barker  pulses.  In  other  words,  when  the 
Barker  pulses  "line  up"  the  entire  data  sequence  is  assumed  to  be 
reasonably  aligned  in  time.  The  Barker  code  was  utilized  so  that  the 
crosscorrelation  function  would  have  a  more  well  defined  peak  compared 
to  some  other  waveform  such  as  a  cw  pulse.  Even  though  much  of  the 
Barker  code  is  distorted  by  the  10  kHz  filters  prior  to  digitization, 
the  Barker  pulse  data  did  exhibit  this  sharp  peak. 

Since  the  filtered  Barker  pulses  are  essentially  bandlimited 
functions,  the  crosscorrelation  of  two  of  them  will  also  be  narrowband 

19 

and  consequently  will  possess  an  envelope  as  a  function  of  time  shift. 

To  minimize  the  error  in  time  shift  required  to  produce  the  maximum 
correlation,  the  envelope  is  considered  rather  than  values  of  the 
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crosscorrelation  function  itself.  Thus  the  alignment  procedure  consists 
of  computing  a  crosscorrelation  envelope  for  several  time  differences 
between  a  given  event  and  a  reference  event  and  shifting  the  data 
accordingly.  A  more  detailed  explanation  of  the  computational  procedure 
is  given  in  the  appendix.  As  a  final  check,  the  phases  of  the  cw  refer¬ 
ence  pulses  for  the  two  data  sequences  are  compared  after  the  data  have 
been  shifted.  In  all  cases  the  resulting  phase  differences  were  less 

than  25°.  The  data  were  further  corrected  for  this  phase  difference 

19 

by  rotating  each  pair  of  quadrature  samples  as  follows: 


X' (na)  =  X(na)  cosAcj)  +  Y (na)  sinA<J> 


and 


Y'(na)  =  Y(na)  cosAtp  -  X(na)  sinA<j>  , 

where  Acp  is  the  phase  difference  between  reference  pulses.  Once  the 
selected  data  segment  had  been  properly  aligned,  it  was  output  to  digital 
tape  for  further  processing.  The  selected  segment  started  at  approximately 
16.0  msec  after  transmission  and  extended  for  6.25  msec. 

In  the  reverberation  data  taken  during  the  1974  experiment,  no 
Barker  pulse  was  available  so  the  crosscorrelation  procedure  was  applied 
to  the  constant  phase  reference  pulse.  Although  the  results  appear 
satisfactory,  the  use  of  a  coded  pulse  should  make  the  procedure  more 
reliable. 

B.  Computation  of  the  Covariance  Matrix  Estimate  and  Data  Segment 
Selection 


The  covariance  matrix  involved  in  the  optimum  detector  was  defined 
previously  by: 


K  =  E{ (N-M) (N-M)1} 
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The  elements  of  this  matrix  are  estimated  from  an  ensemble  average  over 
approximately  100  events.  The  estimated  matrix  is  employed  to  detect  a 
known  signal  added  to  these  same  100  events.  Each  element  of  the  estimated 
matrix  represents  the  covariance  between  the  mth  time  sample  of  the  ith 
element  and  the  nth  time  sample  of  the  jth  element  and  is  computed  by 


VW  -7  Z  “ i 

Y=1  J 


Y) 


|r  VY>J  {7  £  • 


(4.1) 


for  i,j=l,  number  of  elements,  and 
m,n=l,  number  of  samples 

where  n.(t  ,y)  is  the  mth  time  sample  from  the  ith  element  on  the  vth 
l  m 

event  and  r  represents  the  number  of  events  in  the  ensemble. 

This  matrix  may  be  presented  graphically  by  considering  the 
quadrature  representation  as  Frazer  does  in  Ref.  8.  This  representation 
is 


K..(t  ,t  )  =  X..(t  ,t  )  cosa)  [t  -t  ] 
ij  m  n  lj  m  n  o  n  m 


+  Y .  .  (t  , t  )  sinw  f t  — t  ] 
ij  m  n  o  n  m 


where 


X.  .  (t  ,t  ) 
ij  m  n 


■  27  £  [WY)  V£n'Y)  + 


Y.  .  (t 
ij  m 


>'„>  ■  77  z  [wy>  ,'j<h,T)  ■  h(tn'T)]  , 
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and 


wy)  =  xi(ctn-Y)  -  7  i  x;(vy) 

Y=1 


^i(tm*Y)  =  ^(tm*Y)  -  r  E  ^VY> 

Y=1 


with  [x^(t^,y),  f/I(t  *y)1  equal  to,  respectively,  the  in  phase  and 

quadrature  components  of  n.(t  ,y) 

*  1  m 


The  envelope  and  phase  of  the  elements  of  the  covariance  matrix 
are  determined  in  the  usual  manner  from  the  polar  representation 


K..(t  ,t  )  =  E .  .  (t  ,t  )  cosioj  [t  -t  ]  +  4>..(t  ,t  )i 
lj  m  n  ij  m  n  (on  m  ij  m  n  ( 


Thus  the  envelope  is  computed  by 


E.  .(t  ,t  )  =  {x2.(t  ,t  )  +  Y2.(t  ,t  )[ 

IJ  m  n  (  ij  m  n  ij  m  n  ) 


1/2 


(4.2) 


and  the  phase 


S  .  . (t  , t  )  =  tan 
ij  m  n 


-1 


-Y.  .(t  ,t  ) 
n  m  n 


X..(t  ,t  ) 
ij  m  n 


(4.3) 


The  elements  of  the  normalized  covariance  matrix  are  given  by 
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k.  .(t  ,t_)  = 

ij  m  - 


K.  .  (t  ,t  ) 
n  m  n 
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and  the  envelope  of  k..(t  ,t  )  by 

ij  m  n  J 


m* 


t  ) 


n 


E. . (t  ,t  ) 
ij  m  n 


Vi.  .(t  ,t  )  k.  .  (t  ,t  ) 

ii  m  m  jj  n  n 


(4.4) 


Plots  of  normalized  envelope  and  phase  for  the  covariance  matrices  used 
in  this  study  are  shown  in  Figs.  4. 1-4. 4.  While  these  figures  repre¬ 
sent  only  a  slice  (fixed  t  )  of  the  matrix  used  by  the  detector,  the 

functions  should  have  similar  form  for  all  values  of  t  . 

m 

The  particular  time  segment  shown  in  the  covariance  plots  was 
selected  in  the  following  manner.  A  segment  1  msec  long  was  chcsnn  which 
corresponds  to  the  approximate  length  of  the  target  echoes.  Since  one 
intent  of  this  study  is  to  evaluate  the  effect  on  the  detection  problem 
of  considering  the  covariance  between  channels,  it  is  desirable  to  choose 
the  1  msec  segment  in  the  reverberation  data  where  a  significant  covari¬ 
ance  occurs  between  channels.  To  determine  this  point  the  covariance 
function  was  computed  for  zero  time  shift  versus  time  of  observation,  i.e., 

K. . (t  ,t  )  versus  t 
ij  m’  m 

for  each  pair  of  adjacent  channels.  The  envelope  (Eq.  (4.2)),  phase 

(Eq.  (4.3)),  and  normalized  envelope  (Eq.  (4.4))  for  elements  one  and 

two  (K^2 (^ni* are  Presente<^  in  Figs-  4.5  and  4.6  for  the  cw  and  LFM 

data,  respectively.  The  covariance  and  were  similar  to  while 

K, _  and  K_.  indicated  a  much  lower  interchannel  coherence.  For  this  rea- 
45  56 

son  only  elements  1-4  are  considered  in  the  sequel.  The  selected  data 
segment  is  indicated  by  dashed  lines  and  consists  of  320  samples  from 
t  =  18.4  msec  to  t  =  19.4  msec  after  transmission. 


Since  the  effect  of  considering  noise  (reverberation)  on  separate 
channels  is  being  evaluated  here,  it  is  desirable  to  have  the  same  average 
level  of  noise  power  on  each  channel.  If  one  channel  had  a  lower 
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FIGURE  4.4 

PHASES  OF  THE  COVARIANCE  MATRIX 

(1  msec  LFM  REVERBERATION  AT  tm=  18.9  msec  AFTER  TRANSMISSION) 
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FIGURE  4.5 

COVARIANCE  BETWEEN  ELEMENTS  1  AND  2 

(1  msec  cw  REVERBERATION) 
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FIGURE  4.6 

COVARIANCE  BETWEEN  ELEMENTS  1  AND  2 
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reverberation  level  than  the  other  three  channels,  then  this  channel  would 
have  a  more  pronounced  effect  on  the  detection  problem  than  the  other 
three.  Since  this  study  is  attempting  to  evaluate  the  effect  of  inter¬ 
channel  covariance,  rather  than  the  obvious  effect  of  having  little  or  no 
noise  on  one  channel,  a  multiplying  factor  was  applied  to  all  the  returns 
from  each  channel  so  that  the  average  noise  level  is  the  same  as  that  for 
the  channel  with  the  lowest  noise  level.  The  average  used  for  this  pur¬ 
pose  was  the  mean  square  value  of  all  320  samples  from  each  of  the  events 
from  a  particular  channel.  This  process  was  also  performed  on  the  signal 
echoes,  although  in  that  case  there  was  only  one  return  per  channel. 

C.  The  Detection  Problem 


The  detection  problem  considered  here  is  studied  by  adding  a  known 
signal  to  a  series  of  noise  returns.  The  known  signal  consists  of  the 
echo  from  a  styrofoam  sphere  measured  in  a  high  signal-to-noise  environ¬ 
ment.  The  noise  returns  consist  of  reverberation  obtained  by  pointing 
the  transmitter  and  receiver  toward  the  surface  as  described  in  the  pre¬ 
vious  chapter.  To  facilitate  control  of  the  signal-to-noise  ratio  when 
adding  the  signal  to  the  noise  returns,  the  signal  samples  were  adjusted 
so  that  their  mean  square  value  was  equal  to  that  of  the  reverberation 
returns.  Envelope  plots  of  signal-plus-noise  returns  for  channel  one  are 
shown  in  Fig.  4.7  for  both  the  cw  and  LFM  data.  The  dashed  lines  indi¬ 
cate  the  portion  of  the  return  which  is  under  test  by  the  detection  pro¬ 
cess.  In  all  cases  considered  in  the  detector  evaluation,  the  S/N  is  much 
lower  than  the  example  in  Fig.  4.7. 

To  evaluate  the  performance  of  the  detector  developed  in  Section 
II. C. 2,  both  the  predicted  and  the  experimental  ROC  curves  are  presented. 
The  predicted  curve  is  found  in  the  manner  followed  in  Section  II. C. 3 
except  that  under  hypothesis  ,  J/=M+aS,  where  a=10** (SNR/20) .  With  that 
assumption  Eq.  (2 . 12) , becomes 

3=1-  ${<J»_1(l-a)  -  aax)  .  (4.5) 
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where  a  =5  K  as  before.  The  experimental  performance  is  obtained  by 
solving  for  the  vector  of  filter  coefficients  8  (Eq.  2.6)  and  processing 
each  of  the  reverberation  returns.  This  processing  consists  of  computing 
a  test  statistic  X  by  forming  the  scalar  product  of  is  and  W  for  each  mem¬ 
ber  of  the  ensemble  of  reverberation  returns.  A  histogram  is  made  of  the 
resulting  test  statistics 


x  =  nt  b  y  -  l.r 

y  -  y- 


where  y  is  the  event  number  mentioned  earlier  and  T  is  the  total  number 

of  events.  This  histogram  is  formed  by  dividing  the  range  of  X,  i.e., 

(X  -X  .  )  into  20  bins  and  counting  the  number  of  X's  which  fall  with- 
max  min 

in  each  bin.  The  experimental  density  function  under  the  noise-only 
hypothesis  (Hq)  is 


X/H 


(i)  = 


number  of  hits  in  bin  i 


When  the  signal-plus-noise  returns  are  processed  in  a  similar  manner. 


the  density  function  for  X  under  the  signal-present  hypothesis,  p 


a/H. 


(i). 


is  obtained.  Alternatively,  Py/H(i)  may  be  obtained  directly  from 


p. /jj(i)  by  observing  that  under 


X  =  (N+aS)T8  =  NT8  +  aST8 


-*T 

and  hence  p.  .  is  merely  p  .  shifted  rigiit  by  a8  B_.  Points 
A  j  H  ^  A  /  ii 

[  (:t(j )  ,6(j)  J  on  the  ROC  curve  are  then  computed  by 


"(j)  =SPX/H  (i) 

.  .  o 


rj» 


and 


(4.7) 


£(j)  -  P\/H  j  =  1.  20  +  number  of  bins  in  aS^8. 

i>j  1 

Several  comparisons  are  made  in  Che  following  sections  between  the  per¬ 
formance  predicted  theoretically  by  Eq .  (4.5)  and  that  measured  experi¬ 
mentally  by  Eqs.  (4.6)  and  (4.7). 

D.  Analysis  of  Results 


To  evaluate  detector  performance  the  system  of  linear  equations 
indicated  by  Eq .  (2.6)  was  solved  for  several  different  covariance 
matrices.  This  computation  was  performed  on  ARL:UT's  CDC  Cyber  171  com¬ 
puter  using  an  algorithm  titled  LEQT2P  from  International  Mathematical 
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and  Statistical  Libraries.  When  the  size  of  the  system  became  larger 
than  100,  the  linear  system  routine  failed  to  obtain  a  solution.  How¬ 
ever,  for  all  results  presented  in  this  study,  the  solutions  were  accurate 
to  within  10  decimal  digits.  As  an  additional  check  on  the  validity  of 
the  results,  the  variance  of  the  experimental  A's  was  compared  to  that 
computed  by  Eq.  (2.8).  Since  the  covariance  matrix  K  was  computed  from 
the  same  set  of  noise  events  that  were  used  to  test  the  detector,  the 
two  variances  should  be  equal.  This  is  true  in  all  cases  considered 
here . 

1 .  Effect  of  Input  Signal-to-Noise  Ratio 

The  first  results  to  be  presented  are  the  optimum  detector's 
performance  for  various  input  signal-to-noise  ratios  (S/N)  .  Only  one 
channel  is  considered  and  samples  are  taken  according  to  Eqs.  (2.14)- 
(2.16)  with  k=8,  producing  ten  samples  of  each  quadrature  component 
during  a  1  msec  time  interval.  This  represents  sampling  each  quadrature 
component  at  the  Nyquist  rate  since  the  data  are  essentially  limited  to 
a  bandwidth  of  10  kHz.  Figures  4.8  and  4.9  represent  the  predicted  and 
experimental  performance  of  the  optimum  detector  for  the  cw  and  LFM  data, 
respectively.  In  both  cases  the  experimental  and  predicted  values  agree 


PROBABILITY  OF  DETECTION  (0) 


PROBABILITY  OF  DETECTION  (0) 


rather  closely.  Also,  note  that  the  LFM  detector  is  more  powerful  than 
the  cw  detector  for  the  same  S/N.  This  is  generally  the  case  when  the 
transmitted  pulse  lengths  are  equal  and  is  primarily  due  to  the  shorter 
correlation  time  of  the  LFM  reverberation,  i.e.,  the  LFM  covariance 
matrix  is  more  nearly  diagonal. 

2.  Effect  of  Sampling  Frequency 


Next  the  effect  of  sampling  frequency  is  considered.  With 
the  input  S/N  fixed  at  -24  dB,  several  values  of  k  were  utilized  in  the 
detection  problem.  The  results  are  shown  in  Figs.  4.10  and  4.11  for 
k=20,  10,  8,  4,  and  2  for  the  single  channel  cw  and  LFM  data  sets.  This 
range  of  k  results  in  input  vectors  of  lengths  8,  16,  20,  40,  and  80. 

Again  note  the  agreement  between  the  experimental  and  predicted  perform¬ 
ance.  Also,  observe  the  pronounced  improvement  as  the  sampling  frequency 
is  increased.  There  is  a  large  jump  in  improvement  when  k=2,  i.e.,  the 
sampling  frequency  is  40  kHz  per  quadrature  component.  The  close  agree¬ 
ment  between  the  predicted  and  experimental  results  is  indicative  that  the 
phenomenon  is  not  due  to  some  computational  error.  These  results  demon¬ 
strate  that  sampling  frequency  is  not  a  trivial  concern  in  a  study  such 
as  this  one.  Furthermore,  these  results  expose  the  fallacy  of  a  naive 
consideration  of  signal  detection  in  a  bandlimited  environment.  That  is, 
in  both  the  k=4  and  k=2  cases,  the  samples  are  taken  at  a  rate  which 
exceeds  the  Nvquist  frequency.  Thus  one  might  conjecture  that  in  the 
k=4  case  (20  kHz  per  quadrature  component)  one  could  interpolate  with 
the  sampling  expansion  and  thus  reconstruct  the  additional  samples  that 
would  be  present  in  the  k=2  case  (40  kHz  per  quadrature  component) .  Then 
one  might  conclude  that  since  the  Neyman-Pearson  detector  optimally  pro¬ 
cesses  the  data,  the  performance  in  the  k=4  case  should  be  very  close  to 
the  performance  in  the  k=2  case.  However,  as  we  have  seen,  this  conjec¬ 
ture  is  not  true.  The  interpolation  error  associated  with  the  above 
reasoning  is  thus  significant  in  this  context. 
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-  PREDICTED  PERFORMANCE 
EXPERIMENTAL  PERFORMANCE 
O  SAMPLING  PARAMETER  k  =  2 
(40  kHz  QUADRATURE  PAIR 
RATE,  80  TOTAL  SAMPLES) 

O  SAMPLING  PARAMETER  k  =  4 
(20  kHz  QUADRATURE  PAIR 
RATE,  40  TOTAL  SAMPLES) 


FIGURE  4.10 

EFFECT  OF  SAMPLING  PARAMETER  k 
ON  DETECTOR  PERFORMANCE 

(1  msec cw  DATA,  CHANNEL  1.  INPUT  S/N  = -24  dB) 
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3. 


Multiple  Channels 

Most  active  sonars  involve  not  just  one  element  but  an  array 
of  elements.  The  covariance  matrix  associated  with  such  a  multiple  sen¬ 
sor  system  contains  both  spatial  and  temporal  information  about  the  noise 
processes  on  each  element.  The  results  described  in  this  section  relate 
to  the  optimum  detector  utilizing  all  this  covariance  information. 

The  input  to  the  detector  consists  of  four  vectors  v.,  i=l,  4 
as  in  Eq.  (2.1A)  where  each  consists  of  samples  according  to 
Eqs.  (2 . 14)- (2 . 16) .  From  this  point  on  the  detector  is  evaluated  as 
before  in  terms  of  8,  K,  and  S  except  that  and  S  now  represent  groups 
of  vectors  and  K  is  a  matrix  of  submatrices.  In  other  words,  once  B,  K, 
and  are  formed,  they  are  treated  exactly  as  they  were  in  the  single 
channel  case. 

Figures  4.12  and  4.13  demonstrate  the  performance  improvement 
as  the  number  of  channels  is  increased  from  one  to  four.  Here  the 
S/N  =  -24  dB  and  the  sampling  parameter  k=8,  resulting  in  a  sampling  rate 
of  10  kHz  per  quadrature  component  on  each  channel.  As  before  the  experi¬ 
mental  points  agree  closely  with  the  predicted  curves. 

4 .  Suboptimum  Detectors  -  Multiple  Channels 

In  the  following  presentation  comparisons  are  made  between  the 
optimum  detector  and  several  simplified  processors.  The  simplified  pro¬ 
cessors  are  derived  by  making  certain  assumptions  about  the  noise  pro¬ 
cesses  which  result  in  a  simpler  covariance  matrix.  To  predict  the 
performance  of  these  suboptimum  detectors,  the  method  of  Section  II. C. 3 
is  pursued  using  the  assumed  form  of  the  covariance  matrix.  Let 


where  il  is  the  assumed  cover  Lance  matrix  and  8  is  the  array  of  coefficients 
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FIGURE  4.12 

DETECTOR  PERFORMANCE  WITH  NUMBER 
OF  CHANNELS  AS  A  PARAMETER 

(1  msec  cw  DATA.  SAMPLING  PARAMETER  k  =  8, 
INPUT  S/N  =  -24  dB) 
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FIGURE  4.13 

DETECTOR  PERFORMANCE  WITH  NUMBER 
OF  CHANNELS  AS  A  PARAMETER 

(1  msec  LFM  DATA,  SAMPLING  PARAMETER  k  =  8, 
INPUT  S/N  =  -24  dB) 
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for  the  associated  suboptimum  detector.  The  test  statistic  produced  by 

% 

processing  the  data  with  B  is 


% 

X 


with  mean  value 


and  variance 


'V/  rp  'b 

E{X}  =  E{r  }B 


9  'b  9 

a  =  E{ [X-E(X) ]  } 

A 


=  E  j[[£-E{l/}]T!]2j 


'Vr  'v 

=  8iK  8 


As  before  the  matrix  K  represents  the  true  covariance  of  the  noise  pro¬ 
cess.  The  means  under  the  null  and  alternate  hypotheses  are 


*b  %  tA* 

H:m  =E{A}~  MB,  and 
o  o  o  - 


"Vi  «  'n'X'  T»rb 

Hl:ral  =  V*}  =  V-  -  - -} 


=  M18  +  a SLB 


^  T 

=  m  +  aS  8 

o - 


The  performance  is  determined  again  as  in  Section  II. C. 3  to  be 


3  =  1  -  4> 


m  +  rr  4>  ^(l-«) 

o  'S* 
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'p'V 

_l  aS  B 

3  =  1-** 

VB  KB 


(4.8) 


'o  _i 

Note  that,  if  B=K  then  Eq.  (4.8)  simplifies  to  Eq.  (2.12)  as  it 
should.  Also  note  that  the  performance  is  completely  determined  (for  a 
given  S/N)  by 


T% 

.2  A  £  S 
d  =  ,  "• 


VM 


commonly  referred  to  as  the  detectability  index.  The  performance  of 

several  suboptimum  detectors  is  presented  below  by  tabulating  the  values 
2 

of  d  assuming  K=K. 

CASE  I:  Assume  that  the  noise  process  is  wide  sense  stationary 

over  the  1  msec  observation  interval.  In  that  case  the 

covariance  matrix  elements  K..(t  ,t  )  would  only  be  func- 

ij  nr  n  J 

tions  of  t  -t  given  i  and  j.  To  determine  the  effect  of 

n  m  ^ 

this  assumption,  K  is  taken  to  be  K  with  each  element 
replaced  by  the  average  along  the  diagonal  within  its 
submatrix,  i.e. , 


K.  .  (t  -t  )  = 
ij  n  m 


N+n-m 


N  +  n  -  m 


^  Kj j (m-n+k ,k) ,  m  -  n  >  0 


N-n+m 


N  -  n  +  m 


'y  '  (k,n-m+-k)  ,  m  -  n  <  0 


This  operation  produces  a  matrix  which  has  the  same 
general  form  as  the  true  covariance  but  would  appear  to 
have  come  from  a  wide  sense  stationary  noise  process. 
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CASE  II: 


'V 

Assume  the  process  is  independent  in  time.  Here  K  is 
farmed  by  setting  to  zero  all  elements  of  K  except  those 
along  each  submatrix  diagonal  (including  the  main 
diagonal) . 


CASE  III:  Assume  the  process  is  independent  in  space.  In  this 

<\j 

case  K  consists  of  the  autocovariance  submatrices  of  K. 

CASE  IV:  Assume  cases  I-III.  Here  k  is  a  diagonal  matrix  with 

all  diagonal  elements  equal  to  the  average  of  the  diagonal 
elements  of  K.  and  all  other  elements  are  zero.  This  case 
represents  the  assumption  most  typically  made  in  the 
design  of  signal  processors  in  many  modern  sonars.  This 
assumption  leads  to  the  simple  matched  filter  (see  Sec¬ 
tion  II. C. 2)  because  the  filter  coefficients  are  equal 
to  the  signal  samples  within  a  scale  factor  (assuming 
equal  variance  on  all  channels) . 


2 

Table  II  presents  d  for  the  cw  and  LFM  data  sets.  Included 
are  the  optimum  detector  and  CASES  I-IV  for  sampling  parameter  k=8.  When 
the  number  of  channels  in  the  cw  data  is  greater  than  one,  the  linear 
equation  solving  routine  is  unable  to  reach  a  solution  for  the  CASE  I 
problem.  In  all  cases  presented  in  Table  II  the  experimental  performance 
agreed  closely  with  the  performance  predicted  by  Eq .  (4.8). 

The  CASE  I  results  are  especially  significant  in  that  they 
represent  a  detection  performance  of  a  processor  utilizing  a  more  real¬ 
istically  obtainable  covariance  matrix.  A  similar  matrix  might  be 

g 

obtained  from  a  theoretical  model  which  might  accurately  reflect  the 

overall  matrix  structure.  An  actual  covariance  matrix  would  seldom  be 

available.  CASE  I  is  also  significant  because  the  wide  sense  stationary 

23 

assumption  leads  to  a  simpler  processor,  derived  by  Adams  and  Nolte, 
requiring  the  inversion  of  smaller  matrices  with  an  order  equal  to  the 
number  of  channels.  Thus,  the  CASE  I  results  are  more  nearly  representa¬ 
tive  of  those  which  might  be  obtained  using  a  predicted  covariance  matrix 
and  a  simplified  processor  structure. 
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TABLE  II 

2 

COMPARISON  OF  DETECTOR  PERFORMANCE  (d  )  FOR 
OPTIMUM  AND  SIMPLIFIED  PROCESSORS 


r 


5.  Suboptimum  Detectors  -  Single  Channel 

Consideration  of  CASES  I  and  II  for  the  single  channel  problem 
also  yields  interesting  results.  Table  III  presents  the  detectability 
indices  for  CASES  I  and  II,  I  and  II  combined,  and  the  optimum  detector 
for  the  cw  and  LFM  data.  The  sampling  parameter  k  is  varied  from  20  to 
2  corresponding  to  sample  frequencies  of  2.5-40  kHz.  For  both  data 
sets  the  linear  system's  routine  failed  to  reach  a  solution  for  CASE  I 
when  k=2  and  4.  Observe  that  simplified  processors  do  not  exhibit  the 
pronounced  improvement  with  sample  frequency  apparent  in  the  optimum 
detector  results.  The  fact  that  CASE  II  and  the  combination  of  CASES  I 
and  II  are  very  nearly  in  agreement  indicates  that  the  diagonal  elements 
of  the  covariance  matrix  do  not  vary  significantly  over  the  1  msec  time 
interval.  In  other  words,  the  variance  does  not  change  enough  over  the 
observation  interval  to  significantly  affect  the  performance  of  the  sim¬ 
plest  detector  (CASES  I  and  II  combined) . 

6 .  Detection  Probability  versus  Signal-to-Noise  Ratio 

The  degree  of  improvement  of  the  optimum  over  any  of  the 
simplified  processors  may  also  be  observed  by  considering  the  performance 
as  a  function  of  S/N.  For  notational  convenience  Eqs.  (4.8)  and  (4.9) 
are  combined  to  produce 


6  =  1-  4>{<r1(l-a)  -  ad2}  >  (4.10) 

where  again  a=10** (SNR/20) . 

In  Fig.  4.14  6  versus  S/N  is  plotted  with  a=0.1  for  the  4- 
channel  cw  data.  Both  the  optimum  detector  and  the  simple  matched  fil¬ 
ter  (CASE  IV)  results  are  shown  allowing  comparison  of  their  performance. 
At  the  50%  detection  probability  the  optimum  detector  provides  approxi¬ 
mately  a  13  dB  improvement  over  the  simple  matched  filter. 
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COMPARISON  OF  SINGLE  CHANNEL  DETECTOR  PERFORMANCE ( d  )  FOR 
OPTIMUM  AND  SIMPLIFIED  PROCESSORS 


PROBABILITY  OF  DETECTION  (0)  FOR  a  =  0.1 


S/N  -  dB 


FIGURE  4.14 

COMPARISON  OF  OPTIMUM  AND  SUBOPTIMUM  DETECTOR 
PERFORMANCES  FOR  A  0.1  FALSE  ALARM  RATE 

(1  msec  cw  DATA,  4  CHANNELS,  SAMPLING  PARAMETER  k  =  8) 
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Figure  4.15  presents  similar  results  for  the  LFM  data  with  the 
addition  of  the  CASE  I  performance.  Here  the  optimum  detector  provides 
an  approximate  8  dB  improvement  over  the  simple  matched  filter.  The 
improvement  is  less  for  the  LFM  data  than  for  the  cw  because  the  LFM 
noise  process  is  more  nearly  independent  (see  Figs.  4.1  and  4.3). 

The  wide  sense  stationary  (CASE  I)  detector  results  shown  in 
Fig.  4.15  represent  the  degree  of  improvement  that  might  be  obtained  by 
considering  a  reasonably  close  facsimile  of  the  actual  covariance  matrix. 
However,  the  improvement  over  the  simple  matched  filter  is  relatively 
modest  compared  to  the  improvement  the  optimum  system  provides.  The 
CASE  I  improvement  would  most  likely  be  greater  for  the  cw  data;  however, 
the  routine  failed  to  obtain  a  solution  in  that  situation.  In  any  case, 
these  results  again  indicate  the  sensitivity  of  the  detector  to  the 
actual  elements  of  the  covariance  matrix. 
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PROBABILITY  OF  DETECTION  (0)  FOR  a  =  0.1 


OPTIMUM  DETECTOR 


WIDE  SENSE  STATIONARY 
DETECTOR  (CASE  I) 


SIMPLE  MATCHED 
FILTER  (CASE  IV) 


S/N  -  dB 


FIGURE  4.15 

COMPARISON  OF  OPTIMUM  AND  SUBOPTIMUM  DETECTOR 
PERFORMANCES  FOR  A  0.1  FALSE  ALARM  RATE 

(1  msec  LFM  DATA,  4  CHANNELS,  SAMPLING  PARAMETER  k  =  8) 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


In  this  study  several  aspects  of  the  discrete  time  Neyman-Pearson 
optimal  detector  have  been  investigated  using  experimental  acoustic  data. 
The  problem  considered  is  that  of  detecting  a  known  signal  from  an  active 
sonar  operating  in  a  reverberant  environment.  The  noise  background  con¬ 
sists  of  sonar  returns  obtained  by  directing  a  transmitter  and  receiver 
toward  the  wind  driven  surface  of  a  freshwater  lake.  High  signal-to- 
noise  ratio  echoes  from  a  styrofoam  sphere  provide  the  known  signal.  For 
both  the  echoes  and  reverberation,  signals  from  four  elements  of  a  line 
array  were  digitized  for  analysis.  In  this  manner  both  spatial  and  tem¬ 
poral  aspects  of  the  detection  process  were  studied.  The  two  types  of 
sonar  transmissions  employed  are  a  1.0  msec  pulsed  continuous  wave  sig¬ 
nal  with  a  carrier  frequency  of  80  kHz  and  a  1.0  msec  pulsed  linear  fre¬ 
quency  modulated  signal  with  a  10  kHz  sweep  centered  at  80  kHz.  Thus, 
this  investigation  consists  of  evaluating  a  well  known  (but  not  often 
employed)  signal  processing  technique  in  a  particular  application. 

The  detection  problem  addressed  in  this  study  is  a  very  special  case 
of  the  general  detection  problem  encountered  by  designers  of  active  sonar 
systems.  Generally  this  problem  consists  of  detecting  a  contact  of 
unknown  location  given  a  series  of  observations,  i.e.,  multiple  trans¬ 
missions.  Furthermore,  the  contact  echo  structure  exhibits  a  certain 
amount  of  randomness  independent  of  the  noise  background.  In  this  study 
both  the  form  and  location  of  the  signal  are  known  precisely  when  the 
signal  is  present.  For  this  reason  only  a  fixed  time  interval  is  pro¬ 
cessed.  Signal  plus  noise  events  are  created  by  adding  samples  from  the 
echoes  to  samples  of  the  reverberation  returns  at  specified  signal-to- 
noise  ratios.  The  reverberation  is  modeled  as  a  Gaussian  random  process 
with  a  known  covariance  matrix.  For  the  experimental  detector  evaluation, 
this  matrix  is  computed  by  averaging  over  an  ensemble  of  approximately 
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100  reverberation  events.  This  matrix  is  employed  by  a  likelihood  ratio 
detector  to  process  these  same  100  returns. 

While  the  detection  problem  considered  here  is  specific  in  nature, 
it  allows  several  aspects  to  be  studied  under  closely  controlled  circum¬ 
stances.  These  aspects  include:  (1)  performance  versus  input  signal-to- 
noise  ratio,  (2)  effect  of  sampling  frequency,  (3)  type  of  transmission, 
and  (4)  number  of  channels.  In  each  case  comparisons  are  made  between 
predicted  performance  and  that  obtained  experimentally  by  processing  the 
100  reverberation  returns  with  and  without  the  signal  added.  In  all 
cases  the  theoretical  and  experimental  estimates  agreed  rather  closely. 
Comparisons  are  also  made  between  the  optimum  detector's  performance  and 
that  of  several  simplified  processors.  These  simplified  detectors  are 
derived  by  making  certain  assumptions  about  the  covariance  matrix. 
Specifically  the  cases  considered  include  assuming  the  noise  process  is: 
(1)  wide  sense  stationary,  (2)  uncorrelated  in  space  but  correlated  in 
time,  (3)  uncorrelated  in  time  but  correlated  in  space,  and  (4)  un¬ 
correlated  in  both  space  and  time. 

The  results  of  this  study  are  as  follows: 

1.  An  automatic  method  of  alignment  of  multichannel  data  digitized 
from  an  analog  tape  recorder  was  developed.  The  results  indicate  that 
in  all  cases  the  data  were  aligned  by  this  method  to  within  one  sample 
time . 

2.  Experimental  detector  performance  agreed  very  closely  with 
theoretical  predictions  in  all  cases.  Thus,  the  detection  performance 
may  be  pred  d  reliably  from  consideration  of  the  covariance  matrix. 

3.  Gen<_.  illy  the  detection  performance  for  the  LFM  data  was 
approximately  6  dB  better  than  it  was  for  the  cw  data.  In  other  words 
the  LFM  performance  was  the  same  as  the  cw  performance  when  the  cw  data's 
input  sign3l-to-nolse  ratio  was  6  dB  higher  than  that  for  the  LFM  data. 

4.  Defection  performance  for  both  the  cw  and  LFM  data  exhibits  a 
definite  improvement  as  the  sampling  frequency  (and  number  of  samples, 
since  a  fixed  time  interval  is  considered)  is  increased.  While  the 
author  feels  this  result  is  rather  inconclusive,  it  does  indicate  that 
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sampling  frequency  is  a  significant  factor  in  the  detector  performance. 

The  simple  matched  filter  showed  no  improvement  with  increased  sampling 
frequency . 

5.  Using  all  four  channels  the  optimum  detector  provided  the 
following  improvement  over  a  simple  matched  filter: 

(1)  13  dB  for  the  1.0  msec  cw  data,  and 

(2)  8  dB  for  the  1.0  msec  LFM  data. 

6.  As  the  number  of  channels  is  increased,  the  improvement  is  more 
significant  for  the  optimum  detector  than  for  the  simple  matched  filters. 

7.  For  the  LFM  data,  approximating  the  covariance  matrix  in  the 
optimum  detector  by  one  with  similar  form  yields  only  a  3  dB  performance 
improvement  (see  Fig.  4.15)  over  the  simple  matched  filter.  This  improve¬ 
ment  is  modest  compared  to  8  dB  obtained  by  employing  the  exact  covari¬ 
ance  matrix.  Corresponding  results  for  the  cw  data  are  not  available  due 
to  computational  difficulties. 

This  study  might  be  improved  or  expanded  by  considering  the 
following  recommendations. 

1.  A  larger  ensemble  of  reverberation  returns  might  be  obtained 
and  the  covariance  matrix  estimated  from  one  half  of  the  ensemble  could 

e  used  to  process  data  from  the  other  half. 

2.  By  recording  an  ensemble  of  target  echoes,  the  more  general 
problem  of  detecting  a  random  signal  might  be  investigated. 

3.  The  bearing  sensitivity  of  the  detector  could  be  studied  by 
changing  the  apparent  received  angle  of  the  known  signal. 

4.  The  co\  iriance  matrix  might  be  predicted  from  a  theoretical 
model  of  the  reverberation  process.  Thus  the  performance  of  the  detector 
employing  this  matrix  could  be  determined. 

5.  In  an  effort  to  understand  the  unusually  high  performance  in  the 
single  channel,  high  sample  rate  example  (Figs.  4.10  and  4.11),  the  per¬ 
formance  could  be  determined  as  a  function  of  the  number  of  diagonals 
retained  in  the  detector's  assumed  covariance  matrix.  In  other  words 
all  matrix  elements  farther  away  from  the  principal  diagonal  than  n  diag¬ 
onals  could  be  set  to  zero  and  tie  performance  could  be  measured  as  a 
function  of  n. 
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6.  A  solution  for  the  performance  of  the  wide  sense  stationary 
detector  (CASE  I)  could  be  obtained  for  the  cw  data.  This  would  involve 
modifying  the  linear  equation  solving  routine  used  in  this  study. 
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APPENDIX 

COMPUTATIONAL  PROCEDURE  FOR  COMPUTING  THE 
CROSS CORRELATION  ENVELOPE  USED  FOR  ALIGNMENT  OF  DATA  SEQUENCES 


The  following  discussion  describes  the  computational  procedure  used 
to  align  data  from  different  pings  and  channels  and  is  taken  primarily 
from  Ref.  19. 

Consider  two  narrowband  functions, 

f^(t)  =  A1(t)  cos{u>ot+<J)1(t) }  , 

and 

f2(t)  =  A2(t)  cos{u>ot+4>2(t) }  , 

where  A^(t)  and  <j>^(t)  are  the  slowly  varying  (with  respect  to  cosu^t) 
amplitude  and  phase  functions  of  f^(t),  (Jo=2rfo,  and  f^  is  the  center 
frequency.  In  quadrature  form  f^  and  f2  are  given  by 

f1 (t)  =  X  (t)  cosu  t  +  Y.(t)  sinw  t  ,  (la) 

i  i  o  1  o 

and 

f2(t)  =  X2(t)  cosojQt  +  Y2(t)  sinwQt  ,  (lb) 

where 

Xt(t)  =  A^t)  cos».(t)  , 

and 

Yt(t)  =  -A^t)  sin4'i(t) 

The  crosscorrelation  of  f^  and  f2  is  given  by 
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PriEClDlitt  PAC*  IUW-WOT  FUA® 


R12(t) 


f1(t)*f2(t) 


T 


(t)f2(t+T)dt 


(2) 


Inserting  Eqs.  (la)  and  (lb)  in  Eq.  (2)  yields 


R12(t>  ‘  \2M 


COSOJ  T 

o 


+  Y„ 


(t) 


12 


smw  t 
o 


where 


(t)  =  l/2{X1(t)*X2(t)  +  Y1(t)*Y2(t>) 


(3) 


(4a) 


and 


Yr  (t)  =  l/2{X1(t)*Y2(t)  -  X2(t)*Y1(t)}  .  (4b) 

Actually  there  are  additional  terms  in  Eq .  (4)  of  the  form 

(X1(t)  cos2u)ot*X2(t))  which  are  found  to  be  zero  by  the  convolution  theorem. 

To  compute  the  crosscorrelation  between  two  low  pass  functions  such 

as  X1(t)*X  (t),  the  interpolatory  representation  of  this  integral  is  used 

1  1  24 

and  is  given  by 


R'(ma)  =  2a^  W  y]  X^ (na)  X9([nri-n]a)  ,  (5) 

n=l 


where  W  is  the  narrower  bandwidtli  of  the  two  signals  and  a  is  the  sampling 

period,  which  is  less  than  or  equal  to  1/(2W).  Since  the  actual  functions 

* 

X^  and  X2  are  of  finite  duration,  i.e.,  X^(na)=0,  n>N  (where  N  is  the 
number  of  samples  in  the  finite  interval),  Eq .  (5)  reduces  to 

*Strictly  speaking,  this  statement  violates  the  bandlimited  assumption; 
however,  in  practice  Eq .  (6)  is  a  good  approximation. 
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(6) 


2 

l'(ma)  =  2a  W  (na)X2 ( [nri-n ]a) 

n=l 


The  envelope  function  of  is  computed  in  the  usual  manner  by 


E12W  =JfR  (ma)  +  Y  R  (ma)  .  (7) 


Normalizing  the  crosscorrelation  function  is  useful  in  that  it  provides 
a  ready  indication  of  problems  in  the  data,  i.e.,  when  two  similar  func¬ 
tions  are  closely  aligned  the  normalized  value  will  be  near  unity.  This 
function  is  given  by 


Rj  2  t"13) 

pi2 (ma)  yRli(o)  r22(o) 


and  the  envelope  of  p  9  by 


e12(ma) 


E12(ma) 

^Rn(0)  R22(0) 


Combining  Eqs.  (4),  (6),  (7),  and  (8)  yields 


(ma)  1  <1^  X1(na)X2(fnrt-n]a)  Y^(na)Y7(  [m+n]a)  (9) 


41 

+  ^  Xj  (na)Y2([nrtn]a)  -  X2 (na) Y^ (  [m+n]a) 


N  N 

:  <  ^  X21(na)  +  Y^na)  ^  X22(na)  +  Y2?(na) 

n=l  J  Ln=1 
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